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We present a theory of the low-energy excitations of a ferromagnetic metal nanoparticle. In addi- 
tion to the particle-hole excitations, which occur in a paramagnetic metal nanoparticle, we predict 
a branch of excitations involving the magnetization-orientation collective coordinate. Tunneling 
matrix elements are in general sizable for several different collective states associated with the same 
band configuration. We point out that the average change in ground state spin per added electron 
differs from non-interacting quasiparticle expectations, and that the change in the spin-polarization, 
due to Zeeman coupling, is strongly influenced by Coulomb blockade physics. 



PACS Numbers:73.40Gk, 73.20.Dx, 73.23.Hk, 75.60.-d 



The energies of many-particle states in metallic 
nanoparticles and semiconductor quantum dots can be 
measured directly by tracking the dependence of resonant 
tunneling conductance peaks on gate and bias voltages 
0. This technique has been used to study the interac- 
tion physics of quantum dots at weak fields O and in 
the quantum Hall regime [3], and metallic nanoparticles 
in both superconducting i] and normal 0] states. In 
zero-field semiconductor quantum dots at high densities 
and in normal metallic nanoparticles, experiments are 
generally consistent with a model which acknowledges 
electron-electron interactions only in a mean-field elec- 
trostatic term that, because of its long range, gives rise 
to the Coulomb blockade effect |6|. The success of this 
simple interpretation is a consequence of the Fermi- liquid 
character of these interacting electron systems. 

The present work is motivated by recent experimental 
studies of tunneling via discrete energy levels in ferro- 
magnetic Cobalt nanoparticles that find resonance spac- 
ings smaller than predicted in an independent particle 
picture, and a dependence on external field qualitatively 
different from those in paramagnetic metal nanoparti- 
cles. Since bulk ferromagnetic metals have low-energy 
spin-wave excitations in addition to their Fermi liquid 
particle-hole excitations, it is natural, as suggested Q by 
Gueron et al, to seek an explanation in terms of the col- 
lective quantum physics of the magnetic order parameter 
field. In this Letter we address the inter-play of quasipar- 
ticle and collective order parameter excitations in tun- 
neling spectroscopy studies of ferromagnetic nanoparti- 
cles. Our conclusions are based in part on the proper- 
ties of a simple exactly solvable toy model, described in 
the following paragraphs. We conclude that only states 
in which all singly occupied nanoparticle orbitals have 
aligned spins are relevant at low-energies in ferromag- 
netic nanoparticles, and use this observation to derive 
expressions for the tunneling amplitudes of ground and 
excited many-particle collective spin states. Near fields 
where magnetization reversal occurs, many of these low- 
energy spin states have large tunneling amplitudes, and 



contribute significantly to the tunneling spectrum, partly 
explaining the enhanced density of resonances seen in ex- 
periment. 

In ferromagnetic metals, short-range exchange inter- 
actions favor spin-alignment, giving rise to an approx- 
imately rigid spin-splitting of quasiparticle energies in 
the ferromagnetic ground state. Our toy model reflects 
this rigidity by assuming identical exchange constants be- 
tween all pairs of orbitals: 
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In Eq. [I] cj is a nanoparticle orbital energy that incor- 
porates the charge correlation physics neglected in our 
model Hamiltonian and f is the Pauli spin matrix vec- 
tor. The single particle orbitals have an average spac- 
ing inversely proportional to the volume of the nanopar- 
ticle and are expected to exhibit spectral rigidity M. 
The many-particle spectrum of this Hamiltonian follows 
readily from the following observations: i) the total oc- 
cupation of each orbital is a good quantum number ii) 
the interaction term is proportional to the square of the 
total electron spin operator iS tot . The 2 Ns states with 
a given set of iV s singly occupied orbitals have their 
band energy degeneracy lifted by the interaction energy 
— (U/2)Stot(Stot + 1), where the total spin S to t has a 
maximum value N s /2. We show below that only this 
N s + 1-fold degenerate spin-multiplet is relevant to the 
low energy physics of a ferromagnetic nanoparticle. 

We start by considering the ground state of a ferromag- 
netic metal nanoparticle in which majority and minority 
spin quasiparticlcs are occupied up to their Fermi ener- 
gies. We take these to have the values ep a and epi in the 
absence of an external magnetic field at a reference total 
particle number. We are interested in how the nanoparti- 
cle evolves as a function of external field and gate voltage. 
We assume, for simplicity, that the quasiparticle energy 
levels of majority and minority spins are equally spaced 
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near their respective Fermi energies. It follows from the 
considerations of the previous paragraph that the total 
energy, relative to that of the reference state, is 

SE = [e Fa - A/2 + S a /2]SN a + [e Fi + A/2 + Si/2]8Ni 
+ E cb {SN a + 5N l ) 2 /2 

A/2(2S + 1)) + ^-(5, - A/2(25 + 1)) 
A 
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- (V g + 9 fi B B/2)SN a - (V g - gn B B/2)SNi. (2) 

In this energy expression, S a ^ are level spacings, and we 
have added by hand the electrostatic Coulomb blockade 
term, which depends only on the total number of parti- 
cles in the grain. The terms proportional to A in Eq. |^ 
originate from the interaction energy —U Stot(Stot + 1)/2; 
the notational change U — > A/(2So + 1) is motivated 
by the identification, explained below, of A as the spin- 
splitting energy of the quasiparticle bands. Here So = 
(N a o — Nio)/2 is the ground state total spin of the refer- 
ence state and S to t = So + (SN a — SNi)/2. Note that for 
a ferromagnetic particle So will be proportional to vol- 
ume. We have not explicitly indicated the capacitance 
ratio relating the gate voltage and the chemical potential 
of the nanoparticle, and have assumed that only Zceman 
coupling to an external field is important. 

Stability of the reference system ground state requires 
that its energy increases at fixed total particle number 
N both for S to t -> S + 1 (SN a = l,5N z = -1) and 
Stot -> S - 1 (SN a = -l,SNi = +1). From this it 
follows that 

\e Fa - e Fl - A+(S a - Si)/2\ < [S a + S t - 2A/(2S + 1)] . 

(3) 

Since the right hand side of Eq. || ~ N^ 1 , where Na is 
the number of atoms in the grain, it follows that A = 
e Fa — e F i to within a fluctuating mesoscopic correction, 
i.e. that A is the quasiparticle spin-splitting. The energy 
difference between the Stot = (N a — N)/2 states retained 
in our considerations and the Stot < (N a —Ni)/2— 1 states 
we have discarded is A2So/(2So + 1) ~ A, well outside 
the energy range of interest ||. 

It is instructive to minimize the energy by consider- 
ing N a and N as continuous variables, thereby obtaining 
the trend lines around which the mesoscopic ground-state 
spin and charge quantum numbers fluctuate as gate volt- 
age and external field vary. We find that 

C£6N a + C-lSNi = V g + 9f i B B 

C£5N a + CH 1 SN i = V a -gn B B (4) 

where 
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The inverse capacitance matrix C _1 describes the cou- 
pled spin and charge response of the nanoparticle. This 
coupling is responsible for the variation of total particle 
number with external field observed |l^,[ll| by Ono et al. 

Eq. [I| can be solved to obtain the variation of both N a 
and Ni with either V g or B. For example, we find that 
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and that 



dNg/dVg = M u - M ai = Sj - A/2S 
dNi/dV g M aa - M ai S a - A/2S ' 



(7) 
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Level spacings and exchange splittings in metallic fer- 
romagnet nanoparticles can be estimated from spin- 
density functional band structure calculations. For the 
case of Cobalt particles Jl2| S a — 5.55eV /Na, Si — 
l.43eV/N A , 2S = 1.65^, and A/2S = 1.07eV/N A . 
The Coulomb blockade energy is sensitive to the screen- 
ing environment of the nanoparticle; for a grain with 
N A « 1500, Ralph et al. find that E cb > 30meV. 
As expected, the largest elements in the inverse capaci- 
tance matrix are the purely electrostatic Coulomb block- 
ade contributions. Evaluating the right hand side of 
Eq. we find that Stot almost always decreases when 
particles are added. The response of Stot to external 
fields, given by Eq. ^, is smaller than would be expected 
naively, because E cb , the dominant mesoscopic energy 
scale, suppresses charge fluctuations causing the larger 
(majority spin) level spacing to limit this response. Typ- 
ically the relatively large value of E cb suppresses the to- 
tal charge response to both gate voltage and external 
field. For example, for Na = 1500 and E cb = 30meV, 
typical for the particle size studied by Gueron et al., 
we find that SN a = 0.0025dVg + 0.306d(g^ B B) and 
5Ni = 0m02dV g - 0.321d(gfi B B) with energies in meV. 

In the model discussed so far, structure would appear 
in the bias dependence of the tunneling current only on 
the scale of the single-particle level spacing. Following 
Gueron et al, we seek an explanation of experiment by 
invoking spin-orbit coupling and long-range dipole in- 
teractions between the electronic spins. In Cobalt, the 
typical band energy shift due to spin-orbit coupling is 
fill SE ~ lmcV; the dependence of the total band 
energy on spin-moment orientation results from a par- 
tially cancelling sum of spin-orbit energy shifts over all 
singly occupied orbitals. The typical spin-orbit matrix 
element between an individual pair of orbitals in a grain 
is ~ SE/Na, ~ 1/ieV for the grain sizes of interest to 
us. These considerations justify starting from a model 
which neglects spin-orbit coupling of states with differ- 
ent band configurations and uses the bulk magnetocrys- 
talline anisotropy coefficient to describe its effect within 
a single multiplet: 
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FIG. 1. Energy eigenvalues E„ vs. {So,n\S z \So,n) for the 
collective spin Hamiltonian. The curves, offset for clarity, 
correspond to increasing values of a, starting from a = —2 
at the bottom. Here So — 25 and B is oriented at n/4 from 
the easy axis, z. At a = q sw ~ 1.2, the metastable state 
disappears. 



H = n- K(S Z /S ) 2 . 



(8) 



where, for Cobalt [L4j, K ~ 0.073AU meV. Note for 
iVji ~ 1500, if ~ .leV, much larger than either level 
spacing or Coulomb blockade energy scales. 

A theory of tunneling spectroscopy in ferromagnetic 
nanoparticles requires results for tunneling matrix ele- 
ments between many-body eigenstates, found by diago- 
nalizing the collective spin Hamiltonian 



H = £ b and - K{S Z /S ) 2 - gn B B ■ S 



E hand -K/S [Si/S + aB-S 



(9) 
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where K/Sq « 0.05mcV is the natural energy scale 
and a = gusB /(K/ So) ~ 2£?[T] is the (dimensionless) 
strength of the magnetic field. In Fig. [l] we plot the en- 
ergy eigenvalues E n as a function of the expectation value 
of S z for the corresponding eigenstates \So,n) for So = 25 
and B oriented at 8 CKt — tt/4 from the easy axis, z. The 
curves are offset for clarity, starting at the bottom with 
a = —2. For large negative a the ground state spin is po- 
larized along the field direction, and gradually reorients 
toward —z as the field is ramped to zero. At a — a level 
crossing occurs and the ground state state is now polar- 
ized along z. The former ground state, now classically 
metastable and separated from the true ground state by 
a potential barrier, is apparent in the spectrum until the 
classical switching field is reached at asw ~ +1-2. The 
classical down-sweep metastable states appear at nega- 
tive values of a and positive values of {So,n\S z \So,n). 

In order to evaluate the many-particle tunneling ma- 
trix elements, we express these eigenstates in terms of 
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FIG. 2. The tunneling spectral function for the collective 
spin model, at different values of the magnetic field, specified 
by a. The solid and dashed lines refer to the Stot — * So ± 1 
transitions respectively. So = 25 and B is oriented at 7r/4 
from the easy axis z. For < a < q sw ~ 1.2, the system 
is in a metastable state (see Fig. |l]), before the tunneling 
takes place. In (e), a is near q sw , where the metastable state 
disappears and the magnetization reversal occurs. 



microscopic electronic degrees of freedom. For a given 
set of quasiparticle occupations, the microscopic state for 
collective spins oriented in direction il = Cl(9, <fi) is 



i*(^))=ri[ v T(% t ,T 



keD 



(11) 



where D and S are the sets of doubly and singly occu- 
pied orbitals respectively; V-f(O) = cos(0/2) exp(i</>) and 
V x (f2) = sin(6>/2)expH0). 

The matrix element between states with collective spin 
orientations Q,' and f2 for adding an electron with spin- 
orientation f2 cxt to an empty orbital, thereby increasing 
N s and S to t, is 

(*(^ / )|c t (^cxt)|*(^))+ - (Ci'\Ci)^ 2 (n'\n cxt ) 1/2 , (12) 

where (fl'\(l) 1/2 is the inner product of spin-1/2 coherent 
states |15| . Similarly the matrix element for adding an 
electron with spin-orientation — i7 xt to a singly occupied 
orbital, thereby decreasing N s and Stot is 

(*(^')|c t (^ cxt )l*(o)}_ = (n'p)^- 1 ^^)^. (13) 

The tunneling spectral function is defined as 
A±(cj) = J2\(( N ° ± l)/2,n|ct(fi ext )|iV s /2,0)| 2 

n 

xS^lj - (E {Ns±1)/2tn - E Ns/2fi ) , (14) 
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where \(N S ± l)/2,n), £7, 



are collective spin 



eigenstates and eigenvalues of TL in the S t ot = 5o ± 1/2 = 
[N s ± l)/2 manifold respectively. To evaluate the matrix 
elements in Eq. O, we expand these eigenstates first in 
terms of eigenstates of S z and then in terms of spin- 
coherent states jl5| with a definite Cl. Using Eqs. 12, 13|, 
we finally arrive at the simplified expressions 



((N s ± l)/2, n| c t(Aext)|JV s /2, 0} = 2 ^ ^ V CT (fl ext ) 

<r m' m 

x C'(iV s + l) c "(Ay (^A±l m'|^ m ; ± <r\ , (15) 



where c™(N s ) are the expansion coefficients of \N s /2,n); 



are Clebsch-Gordan coefficients 



for adding the angular momenta So = N s /2 and s = 1/2. 

Many-particle tunneling spectra for So — N s /2 = 25 
and cx t — 7r/4 are illustrated in Fig. ^ In these 'up- 
sweep' calculations, we chose \N s /2, 0} to be the true 
ground state for a < and a > aw- In the field interval 
< a < a sw , however, we allowed the system to be ini- 
tially in the metastable stable, discussed in Fig. [p. These 
figures clearly show that the tunneling spectra have a 
regime near the classical switching field, where many ex- 
cited states contribute giving rise to a dense tunneling 
spectrum, like that seen in experiment. 

For a particular ferromagnetic metal particle at a given 
gate voltage, the onset tunneling spectra can be domi- 
nated by an intermediate state with either increased or 
decreased particle number. In either case, the lowest en- 
ergy multiplet could have Stot = <So ± 1/2, although de- 
creases will be more common for particle addition and 
increases for particle removal, as we have discussed. The 
calculations presented here can explain a dense tunneling 
spectra with level spacing ~ K/Sq rather than 6 near 
the classical switching field, but not at general fields. 
The magnetic anisotropy energy landscape of a realistic 
nanoparticle is certain to be more complex than in our 
model calculation. There is clearly a shape anisotropy, 
since the particles used in the experiment are approxi- 
mately hemispherical. A self-consistent calculation flcf l 
on small, irregularly shaped clusters predicts a magnetic 
anisotropy energy that can be much larger than the one 
we have used here. If this were the case, this energy 
scale would no longer be a candidate for explaining the 
experiment. Note however that for ellipsoidal particles, 
the shape anisotropy constant is comparable to the bulk 
constant |14|. On the other hand, magnetostatic interac- 
tions with nearby particles and more subtle surface effects 
Jl7] are probably important but much harder to estimate. 
Finally nonequilibrium effects J^|,0] and spin-orbit cou- 
pling between different band multiplets could also con- 
tribute significantly to the ubiquity of dense tunneling 
spectra observed in experiment. 
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